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We theoretically and experimentally investigate the flow fleld that emerges from a rod-like mi¬ 
crorotor rotating about its center in a non-axisymmetric manner. A simple theoretical model is 
proposed that uses a superposition of two rotlets as a fundamental solution to the Stokes equation. 
The predictions of this model are compared to measurements of the azimuthal and radial microflu¬ 
idic velocity field components that are induced by a rotor composed of fused microscopic spheres. 
The rotor is driven magnetically and the fluid flow is measured with the help of a probe particle 
fixed by an optical tweezer. We find considerable deviations of the mere azimuthal flow pattern 
induced by a single rotating sphere as it has been reported by Di Leonardo et al. [Phys. Rev. Lett. 
96, 134502 (2006)]. Notably, the presence of a radial velocity component that manifests itself by an 
oscillation of the probe particle with twice the rotor frequency is observed. These findings open up 
a way to discuss possible radial transport in microfluidic devices. 
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I. INTRODUCTION 

A key aspect in microfluidics is to understand the the¬ 
oretical fundamentals of flow at low-Reynolds numbers 
and to exploit them directly in applications and devices. 
As the miniaturisation is an ongoing trend, especially lab- 
on-a-chip (LOG) or micro total analysis systems (pTAS) 
provide a fast reaction time, low amount of reagents and 
low waste generation as well as high detection accuracy 
and verifiability [I|. Therefore, the analysis and diag¬ 
nostics of reagents within microfluidic channels require 
switching, driving and mixing of the liquid fluids i- Op¬ 
tical forces, for instance, provide a non-contact trapping 
and manipulation of microobjects or biological cells and 
are sufficient for usage in different microfluidic applica¬ 
tions. The latter are e.g., pumping due to rotation of 
microstructures, sorting cells and intermixing different 
fluids ii- Padgett et al. identify the use of optical 
forces in so-called optical tweezers (OT) as a technology 
with high potential for microfluidic applications in biol¬ 
ogy, medicine or chemistry Q. 

Furthermore, the need of measuring systems within mi¬ 
crofluidic channels with high precision arises. Besides ac¬ 
tuation, systems based on optical forces can be used for 
sensing and characterizing of fluid properties. For exam¬ 
ple, Di Leonardo et al. and Kndner et al. have mea¬ 
sured the flow field in a microfluidic device with OT’s 
In both studies a spherical microrotor has been 
activated by interaction of polarized laser light with a 
birefringent vaterite crystal. The rotation has led to a 
fluid flow and a displacement of an inserted probe par¬ 
ticle. The azimuthal flow field around the microsphere 
has been determined by measuring the shift of the probe 
particle. 


However, no radial component has been detected in 
these studies. Objects in microfluidic systems rotating 
about their center in axisymmetric manner can easily 
produce azimuthal flow, which is also supported by the¬ 
ory 1^. Whether a rotating microstructure can in¬ 
duce a radial flow component at low-Reynolds numbers, 
however, seems to be far less understood. Moreover, 
for axisymmetric rotations, simple symmetry arguments 
suggest that such a radial flow can not occur. In ad¬ 
dition, theoretical works [lol - [l^ on non-axisymmmetric 
rotations appear to be restricted to ellipsoidal bodies 
and prove to be quite challenging from a mathematical 
point of view. The work of Camassa et al. 0 , for in¬ 
stance, treats the problem of an ellipsoid sweeping out 
a double cone in Stokes flow. Here, a rather complex 
spatio-temporal velocity field is obtained and explicitly 
analyzed for certain limit cases (small cone angles, far- 
held behavior etc.). Although a radial flow component is 
predicted, the main focus of this works lies on the advec- 
tion of tracers and the establishment of certain invariants 
for the above mentioned limit cases. Gonsequently, the 
main purpose of this paper is to gain some insight on 
the persistence of such radial flow for the case of non- 
axisymmetric rotations at low-Reynolds numbers. More¬ 
over, the azimuthal velocity held component is investi¬ 
gated and compared to the case of a purely axisymmet¬ 
ric rotation. To this end, a comparatively simple model 
that is based on the superposition of two rotlets MM is 
used in order to approximate the rod-shaped rotor geom¬ 
etry and to make statements about the induced velocity 
helds. 

The mere fact that the rotlet solution as a fundamental 
solution to the Stokes equation for low-Reynolds num¬ 
ber hows basically shares the same velocity prohle as the 


2 


point-vortex solution of the two-dimensional Euler equa¬ 
tion [IMl allows to draw an interesting analogy to a 
recently proposed generalized vortex model M- Here, 
the model coirsists of pairs of like-signed point-vortices 
and it has been shown that the straining of these small- 
scale vortical structures leads to non-vanishing radial ve¬ 
locity components that prove to be of great importance 
for the energy transport in such a two-dimensional flow. 
Although the rotlet and the point-vortex solution apply 
to two very different Reynolds number regimes, this strik¬ 
ing analogy allows to predict a possible radial flow also 
for the Stokes equation. 

In this paper, the predictions of this rotlet model are 
compared to experimental results that allow the detec¬ 
tion of azimuthal and radial flow components influenced 
by microrotors. The experiments were performed using 
a combination of optical and magnetic forces, which en¬ 
ables the manufacture of rod-shaped microrotors and the 
investigation of the flow fields. Recently this method has 
been successfully applied for constructing a micro¬ 
pump for microfluidic channels [T^ . 



Figure 1: Approximation of a dumbbell-shaped body of length 
b and neck length c by a pair of rotlets with strength F fixed 
at ±a. Initially, the body is supposed to be at rest as the sur¬ 
rounding fluid undergoes a rotation with the angular velocity 
u) determining the parameters F and a in terms of the rotor 
geometry b and c. 


II. THE ROTLET MODEL EQUATIONS 


The proposed model captures the main effects of the 
rotors investigated in the experiment. Here, we mainly 
follow the procedure described by Chwang and Wu who 
approximated the velocity field of a rotating dumbbell¬ 
shaped body via a superposition of two rotlets spinning 
about the main axis of symmetry [^ : hence, the result¬ 
ing velocity field only possessed an azimuthal component. 
The main difference in our approach is that we consider 
non-axisymmetric rotations, which turn out to induce an 
additional radial flow component. To this end, we ini¬ 
tially work in a frame of reference where the rotor is at 
rest and assume that the entire fluid is undergoing a ro¬ 
tation about the z axis with the angular velocity to. We 
then consider the velocity field that emerges from two 
rotlets at a and —a 
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as an approximation of the velocity field that is induced 
by the rotor. Note that the last term has been added 
to ensure the correct fluid velocity in the far-held. Fur¬ 
thermore, the rotlet distance a and the rotlet strength T 
have to be considered as free parameters yet to be deter¬ 
mined by the rotor geometry. In order to see how these 
two paramters are hxed, we assume that the dumbbell is 
oriented in ea;-direction as it is depicted in Fig. [TJ As 
a consequence, the velocity held in Eq. o has to sat¬ 
isfy no-slip boundary conditions at x^ = ifoea, and at 
Xc = icBj,, namely u(xf,) = 0 and u(xc) = 0. From 
these boundary conditions, we obtain 


b^ + a^ b 
(&2 - a 2)2 “ 2r 
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In eliminating T from these equations, we obtain a rela¬ 
tionship between the geometrical parameters 



(3) 


which can be solved for a variety of rotor geometries in 
order to obtain the non-dimensional parameter a/h in de¬ 
pendence of c/h. The corresponding solution for a/h de¬ 
termining the rotlet distance from the center is depicted 
in Fig. [21 

At this point, it has to be stressed that the approx¬ 
imation m is not capable of reproducing the no-slip 
boundary conditions over the entire rotor body but is 
only strictly valid at the local points Xb = ±bex and at 
Xc = icBy. However, in the following, we are not so con¬ 
cerned about the exact velocity held near the rotor body 
but rather try to understand the immediate consequences 
of the velocity held in Eq. ([I} on the how pattern in the 
far-held. To this end, we make use of the similarity of 
the rotlet velocity prohle in Eq. o to the point-vortex 
solution of two-dimensional turbulence, which only differ 
by the power law in the denominator (quadratic instead 
of cubic for point-vortices). In this context, it has been 
shown recently that a generalized vortex model [isj ] con¬ 
sisting of pairs of like signed point-vortices leads to a 
radial how component in the far-held. 

Applying these results to the rotlet pair of Eq. (|T|) , a 
multi-pole expansion for |a|/|x| ^ I yields 
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Here, the hrst term corresponds to the velocity held of a 
rotating sphere with radius [2ra;]^/^, which solely leads 













3 


a 


b 



part that does not vanish in average, since (sin^((/3 — w<)) 
is a constant. Thereby the time-dependence of the veloc¬ 
ity field not inherently contained in the Stokes equation 
is due to the inclusion of time-dependent boundary con¬ 
ditions similarly to the case of Taylor’s swimming sheet 
[2lj . The corresponding velocity field determined by the 
radial © and azimuthal ([5]) flow components for increas¬ 
ing a/b are shown in Fig. |3]for t = 0. It can be seen that 
the axisymmetric velocity field of the sphere a/b = 0 gets 
distorted in the x-direction for increasing a/h as the two 
rotlets in Eq. © become more and more remote. 

III. ASSEMBLING, ACTUATION AND 
MEASUREMENTS 


Figure 2: Non-dimensional parameter solution for a/b of Eq. 

in dependence of c/b. In the limit of c/6 —>■ 1, the solution 
of a rotating sphere o = 0 is recovered. 


to an azimuthal component of the velocity field. The 
second term however, can lead to a component of the 
velocity in radial direction, which can be seen in writing 
the multi-pole terms in Eq. o explicitly as 


u(x) = TujBz X 2— — 6^a ■ x -f 15^(a • x)^ — 3— 

l x^^ ’ x^ - 
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The term that allows the velocity field to escape from its 
solely azimuthal shape is the second one, which can be 
seen in projecting onto the radial direction. To this end, 

( cos 

sin ip a I , as 
0 

cos p) sin 9 

well as f = ( sin ip sin 9 | , which yields 
cos 9 


Br ■ [e^ X ea](eo ■ e^) = i sin(2((/j - pa)) sin^ 0. (6) 


Additionally, we change to a frame of reference in which 
the fluid is at rest and the rotor is undergoing rotation, 
i.e. p — pa p — u)t and forget about additional centrifu¬ 
gal forces due to the neglection of inertia in the Stokes 
equation ( p|a; x u| small compared to the viscous force 
A^lV^ul, where p is the density and p is the viscosity of 
the fluid). We obtain 


^2 
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for the radial velocity component and 
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for the azimuthal velocity component. It can be seen that 
the radial component oscillates with twice the rotor fre¬ 
quency, whereas the azimuthal component is a superpo¬ 
sition of a time-independent part and a time-dependent 


To verify the established theoretical model a microro¬ 
tor has to be assembled, actuated and the generated fluid 
flow has to be measured experimentally in a precise man¬ 
ner. The OT is used to assemble chain-shaped magnetic 
microrotors by a simple “pick-and-place” technique . 
Based on a focused laser beam, OT provide non-contact 
trapping, moving and rotation of dielectric microparti¬ 
cles . The force acting on a particle is typically in the 
range of several piconewtons and can be described by 

Fopt = -fcr (9) 

with the trap stiffness k and the displacement of the par¬ 
ticle r from its equilibrium position. A Biotin (B) coated 
transparent microsphere (Si02-Biotin-AR722, 0 = 2.73 
± 0.13 pm, microParticles GmbH) is trapped by OT and 
fused together with Streptavidin (SA) coated magnetic 
microspheres (Si02-MAG-SA-S2548, 0 = 2.61 ± 0.12 
pm, microParticles GmbH) by bringing them into con¬ 
tact. The mutual high affinity of the biomolecules leads 
to a stable binding of the complementary particles and a 
magnetic three-particle rotor can be assembled . The 
embedded transparent microsphere in the middle of the 
rotor provides stable trapping and friction-free rotation 
about the rotational axis of the rotor. For bigger rotors 
magnetic interaction is used to attach magnetic particles 
at both ends of the microstructure (Fig. 0] inset). Al¬ 
though, the outer magnetic particles of the five-particle 
rotor are not connected with the B-SA binding, there 
is no deformation of the rotor up to the rotational fre¬ 
quency of 6.5 s“^. The following experimental investiga¬ 
tion is based on chain-shaped rotors with three and five 
particles. 

A cw laser (Verdi G5, Goherent) with a wavelength of 
A = 532 nm is used and with integration of a spatial light 
modulator (SLM) multiple optical trapping spots can be 
generated (Fig. |3|). Computer generated holograms im¬ 
plemented into the SLM allow simultaneous manipula¬ 
tion of multiple microobjects using a single laser source. 
The SLM based OT’s are commonly referred to as holo¬ 
graphic optical tweezers (HOT) [241 - 1^ . The laser beam 
is expanded to utilize the entire active area of the SLM 
(Pluto-VIS, Holoeye) and the reflected modulated beam 
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Figure 3: Velocity field profiles in the x-y-plane determined by the far-held approximation © and Q for increasing a/fo at 
t = 0 (rotor is oriented in a;-direction). The axisymmetric velocity held of the rotating sphere is recovered for a — 0. As a/b 
approaches its maximnm value determined by the solution of Eq. in Fig. [2] the velocity held becomes more and more 
distorted in the s-direction. 


is projected onto the back aperture of a 100 x 1.4 NA 
microscope objective using two lenses (Li and L 2 ) with 
4-f-configuration. The beam is focused into the sample 
chamber and the microparticles can be trapped in an 
aqueous environment. The sample chamber has a circu¬ 
lar aperture of 8 min diameter, a height of 5 mm and is 
closed on the top and bottom with cover glasses. A high¬ 
speed CMOS camera (USB U1-1225LE, IDS) is used to 
observe the particle position and to analyze it with sub¬ 
pixel accuracy down to 5 nm . The rotation of the as¬ 
sembled microstructures is introduced by a controllable 
external magnetic field, which is generated by four elec¬ 
tromagnetic coils arranged around the sample chamber 
and addressed by a stepper motor driver (TMCM-1110, 
Trinamic). Due to the magnetic torque, the inicrorotor 
follows the external magnetic field and rotates with de¬ 
sired rotational speed. The rotation leads to a flow 
field in the liquid environment, which is measured by in¬ 
troducing a non-magnetic probe particle as illustrated in 


the inset of Fig. SI The probe particle is trapped by a sec¬ 
ond optical spot at a variable distance to the rotor axis. 
The generated fluid flow u leads to a displacement of the 
probe particle which can be calculated using Stroke’s law 
(drag force), 

Fdrag = 6TTr]RpCFU ( 10 ) 

with the fluid viscosity rj (the experiments are performed 
in buffer solution contains approx. 99% water), the probe 
particle radius Rp and the correction factor cp- Due to 
boundary effects of the surfaces the correction coefficient 
has to be implemented according to Faxen’s law [^ . 
which is calculated to be cp = 1.62 for these experiments. 
As the probe particle will remain in the equilibrium po¬ 
sition dehned by the optical force and the drag force, the 
fluid velocity can than be expressed by 


dir ^ Rp cp 
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Figure 4: (color online) Schematic depiction of the holo¬ 
graphic optical tweezer (HOT) setup. The laser beam is mod¬ 
ulated and reflected by the spatial light modulator (SLM) and 
focused with a high-NA microscope objective into the sample 
chamber. The high-speed camera is used to observe and de¬ 
tect the particles in the trapping plane. BE: beam expander, 
L: lens, M: mirror, DM: dichroic mirror. Inset: Microscope 
image of probe particle next to assembled five-particle rotor 
(scale bar is 5 pm). 


The equipartition theorem is used in order to evalu¬ 
ate the trap stiffness k. 




IV. RESULTS 


rotational frequency (1/s) 


Fig. [S]a) shows the displacement of the probe parti¬ 
cle in azimuthal direction by varying the rotational fre¬ 
quency from 1“^ to 4s“^ of a five-particle rotor. The 
rotational magnetic field is applied for approx. 4 s for 
each measurement followed by a pause period to clearly 
identify the distinct regions. Here, the rotor radius is Rr 
= 7.2 pm and the distance of the probe particle to the ro¬ 
tor center is r = 9.8 pm. Note that the measured shift of 
the probe particle is constant for a steady rotation (gray 
shaded areas) and is independent on the instantaneous 
rotor orientation. That indicates that the last term of Eq. 
([8|) is averaged during the measurements. As depicted in 
Fig. [5] b) the azimuthal velocity linearly increases with 
the rotational frequency. The averaged azimuthal veloc¬ 
ities are shown for the five-particle rotor corresponding 
to Fig. [S]a) and a three-particle rotor with Rr = 4.4 pm 
at the distance of r = 7.1pm at different rotational fre¬ 
quencies, exemplary. In evidence, the five-particle rotor 
generates higher fluid flow as the three-particle rotor. 

In Fig. [6] the azimuthal fluid velocity in dependence 
of the distance to the rotor is illustrated for the two dif¬ 
ferent rotor types. The measured points represent the 
average values of five measurements at constant condi¬ 
tions. The graphs reveal that there is a marked decrease 
of the fluid velocity with increasing the distance to the 


Figure 5: (color online) a) The displacement of the probe 
particle increases with higher rotational frequency of the five- 
particle rotor, b) The calculated fluid velocities of five- and 
three-particle rotors show linear behavior. 


rotor. The plotted black lines correspond to fits using 
Eq. ([5]) namely = ar~‘^ + whereas in for a 
rotating sphere only decrease as r~^ has been reported. 

Eq. ([7|) indicates that the radial velocity component 
leads to an oscillating fluid flow with twice the rotor fre¬ 
quency. Fig. Ha) highlights this characteristic by show¬ 
ing the radial displacement of the probe particle relative 
to the period T = 2tt juj. The oscillation is measured by 
a rotor with Rr = 7.2 pm at the distance r = 9.8 pm and 
rotational frequency of 4s“^. The black curve represents 
a sine function fit. The period calculated with discrete 
Fourier transform (DFT) is the half of the period of the 
external magnetic field and thus confirms to the theo¬ 
retical model. The maximum positive displacement of 
the probe particle is reached when the rotor is vertical 
aligned. Vice versa the probe particle experiences its 
maximum negative displacement when the rotor is hori¬ 
zontal aligned. 

At this point it is not possible to calculate the exact 
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Three-particle rotor Five-particle rotor 



distance to rotors center (pm) 

Figure 6: (color online) The azimuthal fluid velocity decreases 
with increasing distance to the rotor. The results are in good 
agreement to established Eq. ((S]) which is illustrated by the 
fits (black lines). 


value of the radial fluid velocity. The constant shift of 
the probe particle in azimuthal direction (cf. Fig. [5] a)) 
leads to an unknown change of optical trapping force. 
Hence, Eq. (HU can not be applied here. However, it 
can be noted that higher rotating frequencies lead to an 
increased amplitude of the displacement (Fig. [7]b)). 


V. CONCLUSION 


this experimental technique provides the measurement of 
fluid flow in microfluidic systems with high precision. 
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In conclusion, we have presented theoretical and ex¬ 
perimental evidence that a non-axisymmetric rotation of 
a rod-like rotor at low-Reynolds number flow leads to a 
quite complex spatio-temporal flow pattern. The validity 
of the established theoretical model is demonstrated by 
means of experimental investigations for azimuthal and 
radial flow helds. In particular, the presence of an oscil¬ 
lating radial velocity coiMonent not generated in purely 
axisymmetric rotations y is a quite intriguing new fea¬ 
ture. It will be a task for future experiments to further 
classify the radial flow component especially with respect 
to its power law behavior with increasing distance from 
the rotor. These findings might be important for mixing 
and driving in microfluidic devices as well as for the the¬ 
ory of swimming at low-Reynolds numbers. Furthermore, 


Figure 7: (color online) a) The oscillating radial displacement 
of the probe particle (green dots) is presented at a rotational 
frequency of 4s“^. The black curve represents a sine function 
fit which reaches its maximum (minimum) when the rotor is 
vertically (horizontally) aligned, respectively, b) The radial 
displacement of the probe particle decreases with increasing 
the distance to the rotor. The dashed lines connect the cor¬ 
responding points for the sake of clarity. 
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